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We propose one way to regularize the fluctuations generated during inflation, whose infrared (IR) corrections 
diverge logarithmically. In the case of a single field inflation model, recently, we proposed one solution to the 
IR divergence problem. There, we introduced new perturbative variables which better mimic actual observable 
fluctuations, and proved the regularity of correlation functions with respect to these variables. In this paper, we 
extend our previous discussions to a multi field inflation model. We show that, as long as we consider the case 
that the non-linear interaction acts for a finite duration, observable fluctuations are free from IR divergences 
in the multi field model, too. In contrast to the single field model, to discuss observables, we need to take 
into account the effects of quantum decoherence which pick up a unique history of the universe from various 
possibilities contained in initial quantum state set naturally in the early stage of the universe. 
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I. INTRODUCTION 



It is widely known that on the computation of the non- 
linear perturbations generated during inflation we encounter 
the divergence coming from the infrared (IR) corrections 
iQ i i i II S B iM [B m m. As the possibility of 
detec ting the non-linear primordial perturbations is increasing 
iTMjjHJl^^ 23, 24, 25, 26, 27, 

IZilMlMinilsllMlMlMlMllj, i becomes more im- 
portant to solve the IR divergence problem for the primordial 
perturbations and to predict their finite amplitude that we ob- 
serve rt3^|39i|4(]t|4r|,|42|l. In our previous work ||42|], we have 
proposed one way to solve this IR divergence problem in the 
single-field inflation model. The key observation is that the 
variables that are commonly used in describing fluctuations 
are influenced by what we cannot observe. This is because we 
can observe only the fluctuations within the region causally 
connected to us. We usually define the fluctuation by the de- 
viation from the background value which is the spatial aver- 
age over the whole universe. However, since we can observe 
only a finite volume of the universe, the fluctuations evaluated 
in such a way are inevitably influenced by the information 
contained in the unobservable region. In general, the devia- 
tion from the global average is much larger than the deviation 
from the local average, which leads to the over-estimation of 
the fluctuations due to the contribution from long wavelength 
fluctuations. In addition to that, to discuss the so-called ob- 
servable quantities in the framework of the standard cosmo- 
logical perturbation (even though people call it gauge invari- 
ant formulation), in practice it is necessary to fix the gauge, 
say, the flat gauge. As long as the gauge is determined by 
solving elliptic-type equations on each constant time slice, the 
gauge choice is inevitably affected by the information in the 
region causally disconnected to us. Gauge dependence of the 
perturbation variables itself is not a problem since the "true 
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observables" such as the statistical property of the sky map 
of the temperature fluctuation of the cosmic microwave back- 
ground are not affected by the gauge choice. However, if di- 
vergences appear in the quantities computed at the intermedi- 
ate steps such as n-point functions of the perturbation fields, 
it becomes almost impossible to extract information about the 
"true observables" from them. 

To shut off the influence from the unobservable region of 
the universe, we focused on the presence of residual gauge de- 
grees of freedom in flat gauge. To remove the harmful part in 
the residual gauge degrees of freedom, we imposed a further 
gauge condition which insists the local average of the infla- 
ton to vanish. Then, the fluctuation is not influenced by the 
information from the causally disconnected region. We gave 
a proof that IR divergences are absent in this new scheme. In 
this paper, we extend our discussion to the multi-field models. 
Since the adjustment of the average value is possible only for 
one field, even if we adopt the local gauge, the fluctuations 
of the other components of the scalar fields (p^ are still influ- 
enced by the causally disconnected region. Here we denote 
the index for the D-component scalar fields by / = 1, • • • , D. 
Reflecting this fact, when plural fields are associated with the 
IR divergences, the prescription presented in our previous pa- 
per |42] is not enough to regularize IR divergences. To remove 
the influence from what we cannot observe, we introduce new 
perturbation variables from which we can compute the "true 
observables". As the new variables, we consider the local per- 
turbations of the fields defined by the deviation from their lo- 
cal average values. Then, we prove the regularity of n-point 
functions of these new variables. 

We should note that n-point functions for the local pertur- 
bations are still influenced by what we cannot observe. This 
is due to the difference between the quantum state of the uni- 
verse which we set as a natural initial condition and the one 
which we observe in our real universe. The natural wave 
function of the inflationary universe is not peaked at a spe- 
cific point in the space of the local average values of (p-' . At 
the observation time, this "natural" state of the universe can 
be decomposed into a superposition of wave packets which 
have a peak at a certain point. As the universe evolves, con- 
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stituent wave packets lose correlation to each other Through 
this so-called decoherence process, the coherent superposition 
of the wave packets starts to behave as a statistical ensemble of 
many different worlds, where each world means the universe 
described by a decohered wave packet ll43i 0, IHtl- Our ob- 
served world is just a representative one expressed by a wave 
packet randomly chosen from the various possibilities. Once 
we select one wave packet after the decoherence occurs, the 
evolution of our world will not be affected by the other par- 
allel worlds. However, the initial quantum state does include 
the contributions from all the wave packets. This implies that 
a naive computation of n-point functions is contaminated by 
the other worlds uncorrected to ours. In this paper, taking into 
account the decoherence of the quantum state of the universe, 
we propose a way to define "observables" and show that they 
are actually finite without suffering from IR divergence. 

However, to be honest, the "observables" that we introduce 
do not correspond exactly to what we measure in the actual 
observations. Not to mislead the reader, we should stress here 
that our definition of "observables" does not respect the aspect 
explained in the preceding paragraph in a completely satisfac- 
tory manner They are not the expectation values for a single 
decohered wave packet. What we define is not completely free 
from the contamination of the other worlds decohered from 
ours. However, since we define our "observables" so as to 
over-estimate the amplitude of fluctuations, the proof of their 
finiteness ensures the finiteness of what we really observe. Re- 
cently the stochastic approach |45, 46] has been employed in 
order to solve the IR divergence problem ll39i I40ll41[| . This is 
in harmony with our claim. In stochastic approach, we assume 
that the modes that exceed a certain length scale are automat- 
ically decohered [41 gl gl gl gl |H IH E?]. Namely, 
long-wavelength fluctuations are treated as the statistical vari- 
ance. Since our unique "world" is one realization in this statis- 
tical ensemble, there is no contribution to the observed quan- 
tum fluctuations from long-wavelength modes by definition. 
This is a practical way to take into account the quantum de- 
coherence in the inflationary universe, but this scheme cannot 
completely remove the artificial IR cut-off scale from the dis- 
cussion. Moreover, it is hard to deny the spiteful criticism that 
the reason why the problem of IR divergence does not appear 
in the stochastic approach might be simply because quantum 
fluctuations in the IR limit are neglected by hand. In contrast, 
in our approach, to avoid under-estimating the amplitude of 
fluctuations, we accepted small contamination from the other 
parallel worlds, which will make the amplitude of fluctuations 
larger Namely, sacrificing the accuracy of the estimate of the 
amplitude of fluctuations, we choose to show the IR regularity 
of the observables by over-estimating the amplitude of fluctu- 
ations. 

This paper is organised as follows. In Sec. |ll] we give 
the set up of our problems. Following it, we propose a pre- 
scription to define "observables". The basic idea of our pro- 
posal that assures the IR regularity is stated in this section. In 
Sec. |lll]we explain the details of the proof of IR regularity In 
Sec. |IV] we summarize our results. On the computation of the 
non-linear corrections, the both integrations over the temporal 
and spatial coordinates can make IR corrections singular. In 



this paper, maintaining the initial time ti at a finite past, we 
restrict ourselves to the time evolution of perturbations for a 
finite period of time during inflation. In Sec. IIVI we add the 
comments on the case in which we send ti to distant past. 



II. A SOLUTION TO IR PROBLEM 



A. Setup of the problem 

We first define the setup of the problem that we study in this 
paper We consider the multi-component inflation model with 
the conventional kinetic term. The total action is given by 

(2.1) 

where Afpi is the Planck mass. We denote the field-space met- 
ric by Qij. For simplicity, we assume Qu is a constant ma- 
trix. We perform the following change of variables 



to factorize Mpj from the action as 



(2.2) 



M r 

j V~9[R~ .g'^^/.M'/'f. - '^Vmd^x , (2.3) 

where 4>i = Gij(f'^ ■ Since the Planck mass is completely fac- 
tored out, the equations of motion do not depend on it. The 
Planck mass appears only in the amplitude of quantum fluc- 
tuation. Namely, the typical amplitude of fluctuation of $^ is 
H, and hence that of 4>' is H /Mp\. 

In order to discuss the nonlinearity, it is convenient to use 
the ADM formalism iUgI [l9i I20I1 . where the line element is 
expressed in terms of the lapse function N, the shift vector 
N"-, and the spatial metric hab'- 

ds^ = -N'^dt^ + habidx" + N°'dt){dx'' + N'^dt) . (2.4) 

Substituting this metric form, we can denote the action as 

A'P f r 1 
5* = — Vh N^^^R-^ 2NV{<j)) + —{EabE"^ - E^) 



N 



1 

N 



Nh''''da<i>idb(t)^ d^x , (2.5) 



where 



Eab = ^{hab - DaNb " DbNa] , (2.6) 

E = h'^^Eab , (2.7) 

and D is the covariant differentiation associated with hab- A 
dot " ' " represents a differentiation with respect the time co- 
ordinate. In the ADM formalism, we can obtain the constraint 
equations easily by varying the action with respect to N and 



3 



7V°, which play the role of Lagrange multipliers. We obtain 
the Hamiltonian constraint equation and the momentum con- 
straint equations as 

~2V - N-'^{E''^Eab - E^) 

- h''''da^'db4>i = , (2.8) 
Da[N-\E\ - S'^^^E)] ~ N-^dbM^' - N'^da^') = . 

(2.9) 

Hereafter, neglecting the vector perturbation, we denote the 
shift vector as Na ~ daX- In this paper we work in the flat 
gauge, defined by 



ah 



(2.10) 



where a = e'' is the background scale factor Here we have 
also neglected the tensor perturbation, focusing only on the 
scalar perturbation, in which the IR divergence of our interest 
arises imH. 

In this gauge, using N, x and the fluctuation of the scalar 
fields Lp^ , the total action is written as 



S 



Mil I ^ ^ 
I dtd^xe^P 



n=0 



■■if 



+ N-^{-6p'' + Ap Ax 

-A^Va(^/V>^] , (2.11) 
and two constraint equations are 

2iv2^lF,,...,„(0)/i---/"-6p2 

n=0 

+ ApAx + {V^V^xVaVfeX - (Ax)'} 

+ (07 + V'^xVaV'/)(0' + 0' - VVV,(^^) 

+ N^Va^lV'^f' = 0, 

(2.12) 

{VaN){2pS% + (V«VbX - ^Ax)} 

-{Vbipi)N {4,' + - V^xVa/) - , 

(2.13) 



where 



represents the three dimensional partial differentiation with 
respect to the proper length coordinates X = ePx and 



Spatial indices, a, 5, • • • , are raised by 5"^^. We use this nota- 
tion, which respects the proper distance, because it eliminates 
all the complicated scale factor dependences from the action. 
We define the derivatives of the potential as 



,(0) 



a"T/(<?!)) 



(2.14) 



The background quantities p and satisfy the following equa- 
tions : 



4>' + 3p</.^ + = , 



(2.15) 

(2.16) 
(2.17) 



Expanding the variables as 



1 



N = l+5Ni + -SN2 
1 

X = Xi + H , 

^ = ^1 + H > 



we find that the first order constraint equations are 

Viifii + 2V6Ni + 2pA^xi , +<Pi'pi = , (2.18) 
V, (2p (JA^i - (^/V'f) = . (2.19) 

The first order perturbation ip{ is identified with the field per- 
turbation in the interaction picture. Taking the deviation of 
the action with respect to (pj, we can derive the equation of 
motion for ipi, which includes the Lagrange multipliers SN 
and X- For example, from the third order action, we can de- 
rive the equation of motion with quadratic interaction terms as 
follows. 



+ ipip^ - A/ + Vjif' - 0^Ax + 5NV^ - 3/5 4>^5N ~ dt{6N^^) 

1^5^^^^^ - V„(^^ - 0^(57V)V'^X - {v' - 'i>'SN)A\ ~ p\7-x^a^' ~ ^^(V'^xVa/) 



(2.20) 
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Solving the constraint equations at each order, we can express 
the lapse function and the shift vector as functions of ip^ at 
lower order. Substituting these expressions into the equa- 
tion of motion for (f^, which is up to third order given by 
Eq. (12.20b . the equation is written solely in terms of the dy- 
namical degree of freedom, ip^ . 

B. Tree-shaped graphs 

In this subsection, as a preparation for computing n-point 
functions of (pi{x), we consider an expansion of the Heisen- 
berg field ipi (x) in terms of the interaction picture field, us- 
ing the retarded Green function GRj(a;, x'), which is causal. 
Since the retarded Green function Gyij{x, x') has a finite non- 
vanishing support for fixed x and t' , its three dimensional 
Fourier transform with respect to x' becomes regular in the 
IR limit. 

Let us denote the equation of motion for Lp^ schematically 

by 

C'jp' = ~T'[p], (2.21) 

where £ j is a second order differential operator correspond- 
ing to the linearized equation for (Eq. (I2.32l l) and 
stands for all the nonlinear interaction terms. Using the re- 
tarded Green function 6*^^(2;, x') that satisfies 

C-'jGkk{x, x') = -a-H\x - x')5'k , (2.22) 
we can solve Eq. ( |2.21| i formally as 

ip\x) = ip{{x) + j dy G^^j{x,x')a^{t')r^[ip\{x'){2.23) 

where the first order perturbation Lp{ satisfies 

C'jpi{x) = Q. (2.24) 

Here the factor originates from the background value of 
yj—g. Substituting the expression ( 12.231 1 for p'' {x) iteratively 
into Ti[p\ on its r.h.s., we obtain the Heisenberg field (p^ (x) 
expanded in terms of p{ {x) to an arbitrary high order using 
the retarded Green function G^^j{x, x'). As we have shown 
in ll42ll . to expand p^ , a diagrammatic illustration will be use- 
ful. The Heisenberg field can be expressed by a summa- 
tion of tree-shaped graphs in which all the retarded Green 
functions GYijix^x') are followed by two or more p{{x') 
or GYij{x' , x") with some integro-differential operators and 
all the interaction picture fields p{{x) are located at the right 
most ends of the graphs. 

When we compute the expectation value for n-point func- 
tions of the Heisenberg field, the interaction picture fields 
ipi in the tree-shaped graphs are contracted with each other 
to make pairs, which are replaced with Wightman functions, 
G^^ix.x') EE {p{{x)pi{x')) orGlJ{x,x'){= Gi^{x',x)). 
These propagators are IR singular (oc 1/fc^), which is the 
possible origin of IR divergences in momentum integrations. 
While, the retarded Green function 

Gr^(x, x') = ie{t - t')Ml,{G+'j{x, x') - GJj{x, x')} , 

(2.25) 

is regular in the IR limit. 



C. Iteration scheme and local gauge conditions 

In our previous paper |42], we have shown that the flat 
gauge still has residual gauge degrees of freedom. For in- 
stance, we can introduce an arbitrary function fn{t) to the 
n-th order lapse function and the shift vector as 

bp 

This gauge degree of freedom corresponds to the scale trans- 
formation. As is mentioned in § U our final goal is to define 
finite observable quantities in place of the naively divergent 
quantum correlation functions. For this purpose, we need to 
define gauge invariant variables without the information con- 
tained in the region far outside O. Then, we have to fix the 
residual gauge only using the information near the observable 
region O. 

In the multi-field model, it is convenient to decompose the 
perturbation into the adiabatic one, which is tangential to the 
background trajectory, and the entropy one, which is orthog- 
onal to the background trajectory ll53ll . Using the residual 
gauge degrees of freedom, we fix the homogeneous mode in 
the direction of the background trajectory e/ = 2 
as 

Wt eip\t) = j^JSx Wt{x) ei^\t, x) = 0, (2.26) 

where Wt {x) is a window function, which is unity in the finite 
region Ot = O n with a rapidly vanishing halo in the 
surrounding region, where means a t =const. hypersurface 
corresponding to the time t. For definiteness, we introduce 
D Otj and define O' as the causal past of O'^^. We 
require Wt {x) to vanish in the region outside O'. In addition, 
Wt{x) is supposed to be a sufficiently smooth function so that 
an artificial UV contribution is not induced by a sharp cutoff. 
Lt, an approximate radius of the region Ot, is defined such 
that the normalization condition 

Wtl = l, 

is satisfied. 

We associated " ~ " with the variables in the particular gauge 
satisfying Eq. ( |2.26l l, in order to clearly distinguish them from 
the variables for which the additional gauge condition is not 
imposed. The difference between the variables with and with- 
out " ~ " is only in the boundary conditions. Hence, they obey 
the same differential equations, ( |2.19I )-( I2.20I ). 

In order to fix the arbitrary functions /„(i) {n = 
1, 2, 3, 4, • • • ) so as to satisfy the gauge condition ( 12.261 ). we 
need to obtain a formal solution for p. The higher order lapse 
functions are determined by solving the momentum constraint 
given in the form 

Va (^SN,, - Yp^iv'^ = si") , (n = 1, 2, 3, • • • ),(2.27) 

where the r.h.s. is the n-th order nonlinear term expressed in 
terms of the lower order lapse functions, shift vectors, and p. 
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As we neglect the vector perturbation, we consider only the 
scalar part of these equations, i.e. its divergence, which is 
formally solved as 



6Nn = 6Nn + fn 



(2.28) 



with 



2p 



We define the operation A ^ by 



in 



Wt{e-PY)(fY 
\X-Y\ 



F{t,e-PY) ,(2.29) 



so that it is completely determined by the local information 
in the neighborhood of Ot- Similarly, the higher order shift 
vectors satisfy the Hamiltonian constraint in the form 



p 

V 



--(/„ + A-iV"S(")j +C„, 

where C„ on the rh.s. is a function expressed in terms of the 
lower order lapse functions, shift vector and <^^. A formal 
solution for Xn is given by 



Xn — Xr. 



6p 



■ fn J 



(2.30) 



with 



Xn 



A- 



4>i - -^dt (pifi) 



Substituting the expressions for the lapse function ( 12.28b and 
the shift vector (12.301 1 into the equation of motion for trun- 
cated at the n-th order, we obtain an equation 



p 



2V' ] fn = ~Wt{x)Ti, (2.31) 



where, for later convenience, we have inserted a window func- 
tion Wt{x) on the r.h.s., which does not alter the evolution in 
O. The explicit form of C'j is given by 



C'j^(d^+3pdt-A)S'j+{v!, 



(2.32) 



with 



and on the r.h.s. of Eq. ( |2.31| l represents all the n-th order 
nonlinear terms expressed in terms of lower order variables. 



The equation for the inhomogeneous part of ipl^ is obtained 
by acting the operator 1 — Wt on Eq. ( |2.31t as 

(1 - Wt)Cijf>i = -(1 - Wt)Wtix)Ti,nm . (2.33) 
Then, we find that (pj^ (x) is obtained by 



where W/, = 5j — e^ejWt and ipf^ satisfies 



(2.34) 



(2.35) 



^iri(^) ^ homogeneous field perpendicular to e^. The so- 
lution ( |2.34| i satisfies the gauge condition Wt ei(pn — 0. 

The remaining unknowns are /„ {t) and B^^^ (t) which has 
D — \ components. These D unknown components are de- 
termined by the homogeneous part of the equations of motion 
obtained by substituting ( |2.34t into Eq. ( 12.311 ), 



^ jB±n - fn- 



fn = C\je^eKWt^n- 

(2.36) 



D. Projection to one decohered wave packet 

When plural fields have scale invariant or even redder spec- 
trum, the entropy perturbation can give divergences. However, 
in this subsection we show that a naive computation of the cor- 
relation functions does not give the correlation functions that 
we actually observe. 

When there is no isocurvature mode related to IR diver- 
gence, making use of the gauge degree of freedom, we can 
arrange that the adiabatic perturbation variable A = ei(p^ to 
be the deviation from the local average value. In contrast, 
there is no such gauge degree of freedom for the isocurvature 
perturbation 



e ejip 



Hence, we have to use the isocurvature perturbation variables 
defined by the deviation from the average values on a whole 
time slice, which contains information of the causally discon- 
nected region. As observable isocurvature perturbation, we 
introduce the local fields. 



S' [x) = S' {x) - WtS\x) . 



(2.37) 



However, even if we restrict our attention to the local quantity 
Si{x) on the final surface, the variables Si{x) which contains 
the information outside the causal region appear in describing 
time evolution of the field. Although in our previous work 
we have stressed that the locality of the observables is the 
key issue in order to assure the IR regularity, the locality is 
inevitably violated under the presence of IR divergence origi- 
nating from isocurvature perturbation. 

Here we need to raise another key issue, i.e. the quantum 
decoherence. The primordial perturbations are expected to 
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decohere through the cosmic expansion and/or through vari- 
ous interactions ll43i l44ll45ll . This decoherence process trans- 
mutes the quantum fluctuations at a long wa velength to sta- 
tistical variances S gll Ei Si Isd lSlL El Bltl. At the 



initial time when the wavelength of relevant modes is short, 
the adiabatic vacuum state will be a natural vacuum state. 
However, the adiabatic vacuum state is not a wave packet 
sharply peaked around a specific value of the homogeneous 
part of the scalar field Wt(p^ . Instead, it is infinitely broad 
and can be interpreted as a coherent superposition of peaked 
wave packets. (Detailed explanation will be given in § IB II 
below.) In the early stage of inflation these wave-packets cor- 
relate to each other, but the quantum coherence is gradually 
lost in the course of time evolution. Thus, at the observation 
time it = tf) the coherence will remain only between adja- 
cent overlapping wave packets. Our observed world is corre- 
sponding to one decohered wave packet picked up from this 
superposition. For the later time evolution of our world, we 
can completely neglect the other wave packets whose peak is 
located very far from ours in the space of isocurvature com- 
ponents of the local average values of fields = {Wtfip'')±- 
Hence, keeping the superposition of all wave packets as the 
wave function of the universe gives rather misleading results, 
i.e. huge over-estimates of quantum fluctuations. We should 
therefore remove the contributions from the other worlds. 

It is standard to discuss the decoherence process by coarse- 
graining some degrees of freedom in the quantum interacting 
system, by which the reduced density matrix evolves from its 
initial pure state to a mixed state. This process is interpreted as 
the transition from the initial coherent superposition of many 
different worlds to the final statistical ensemble of them. In 
ll39il40ll4lll . the decoherence process of the long wavelength 
modes is treated by using the stochastic approach to infla- 
tion, in which all the quantum fluctuations of long wavelength 
modes are assumed to turn into the variance of classical en- 
semble at each time step. This assumption of complete clas- 
sicalization can be justified to some extent by coarse-graining 
the short wavelength modes. On physical ground, we believe 
that this approach gives a good approximate description of the 
dynamics of inflation. [15 811 However, here we take a different 
approach because in the stochastic approach, by assumption, 
the quantum fluctuations of long wavelength modes, which we 
focus on in the present paper, cannot enter into the quantum 
loop corrections from the beginning. In this regard we think 
that stochastic approach is not much more satisfactory com- 
pared with a naive prescription of introducing a cutoff length 
scale by hand. 

The accurate evaluation of what we really observe requires 
to elucidate the decoherence process of the primordial per- 
turbations, which is a long-lasting and unsettled issue. Here, 
instead, we aim at proving the IR regularity of our "observ- 
ables" independently of the details of the decoherence pro- 
cess, which is the heart of this paper Although it is difficult 
to understand how the decoherence process proceeds until the 
observation time tf, it is natural to expect that the wave func- 
tion of the universe has been already decohered at t = tf 
to a large extent. The observation picks up one world from 
the superposition of many decohered worlds. The wave func- 



tion corresponding to each decohered world wifl have a rather 
sharp peak in the coordinate space of 



S''it) = e'}WtS\x), (a = 2, 3, 



.D) 



(2.38) 



where {e"} with e} = e/ is a set of orthonormal bases in field 
space. Hence, we insert a projection operator V^a} which 
restricts the values of 5"(t/) to a small range near — 
a" without making any significant effect on each decohered 
wave packet. Then, the insertion of V is expected to reduce 
the contamination from the other worlds significantly. 

For simplicity, we choose the Gaussian projection operator 



V{a} = exp 



CjiS^itf)-a^){Sf'{tf)-a'') 



(2.39) 



where a" are _D — 1 real C-numbers ' . The dispersion should 
be sufficiently large compared with the width of one deco- 
hered wave packet to guarantee that the evaluated amplitude of 
fluctuations is always larger than that for a single wave packet. 
Inserting an identity 



1 



VdetC 



n 



(2.40) 



we can expand the n-point function of the variables whose lo- 
cal average values are subtracted. {A,S'} are schematically 
denoted by O. We can expand the n-point function of 0{x) 
for the adiabatic vacuum I ) „ as 



a{0\d{tf,Xi)d{tf,X2) ■ ■■ditf,Xn)\0)a 
' D 

n 



1 



VdctC 



.a=2 



da" 



/2n 



Xa{0\V{a}O{tf,Xi)---O{tf,Xn)\0)a ■ (2.41) 

Then, we regard 

(V d{tf,xi)ditf,x2) ■ ■ ■ d{tf,x„)) 

_ a{Q\rditf,Xl)d{tf,X2)---d{tf,Xn)\0)a 



a{0\r\0)a 



(2.42) 



in the integrand on the right hand side of Eq. ( 12.411 ) as the ob- 
servable n-point function of Os after the selection of a single 
decohered world. Here we set a" = and denote PjaQ^o} 
by V. Setting a" = does not lose generality because the 
classical average values of isocurvature perturbation can 
be changed by choosing the background trajectory. We will 
prove the IR regularity of this n-point function in the succeed- 
ing section. 



' Here we used the word "projection operator", but this operator does not 
satisfy the relation T-'^^j = "P^^} expected from its name. However, this 
kind of property is unnecessary for our present discussion. 
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Now the question is how to determine the width of the 
projection operator, a. (Here we are assuming that C"'^ ~ 
a^S"^ .) On one hand, a must be chosen large enough to ex- 
ceed the width of a decohered wave packet. Naively there 
is a minimum size of the wave packet since a very narrow 
wave packet cannot maintain its width for a long period of 
time. Later, we find that the minimum size of a wave packet 
that we can choose is determined by the typical amplitude 
of quantum fluctuations generated during inflation, which is 
characterized by the Hubble scale for a nearly massless scalar 
field. This amplitude is H/Mp\ in terms of the fluctuation of 
0{x). Therefore we need to set a to be larger than H/Mp\. 
On the other hand, in order to suppress the influence from 
other wave packets, a should not be very large. Later, we 
find that the condition that the higher order contributions 
are more suppressed requires cr to be much less than unity. 
These two conditions are compatible by choosing a to satisfy 
H/Mpi < o- < 1. 

Due to the inaccurate evaluation of the decohered wave 
packet, the effect of insertion of the Gaussian projection is 
not equivalent to selecting our world through the actual de- 
coherence process. Hence, we cannot claim that the n-point 
function given by Eq. (I2.42l i is the true observable n-point 
function. However, the former amplitude is larger than the 
latter one. Thus, if the n-point function given by Eq. ( |2.42t 
is proved to be finite, we can conclude that the n-point func- 
tion of 0{x) evaluated for the actual decohered wave packet 
is also finite. 

In the above discussion we assumed that the average val- 
ues S°'{t) of all entropy modes has accomplished decoher- 
ence process successfully before we measure n-point func- 
tions of 0{x) at i = tf. Here we want to stress that whether 
the superposition of decohered wave packets come to be sta- 
tistical ensemble or not has nothing to do with whether the 
mode is measurable for us or not. Let's consider a hidden 
variable x which interacts extremely weakly with our visi- 
ble sector Even in that case, if x represents an average of 
a field over a large volume, the quantum coherence between 
two wave packets |1) and |2) peaked at largely different val- 
ues of X will be lost (at least after integrating out the other 
degrees of freedom in the hidden sector). Assuming that x 
takes the two discrete state |1) and |2) with an equal weight 
for simplicity, the evolved density matrix will be schemati- 
cally written as p = (|l)(l|pi + \2){2\p2)/2, after integrating 
out the other degrees of freedom in the hidden sector Here pi 
and p2 are the density matrices of our visible sector (If the in- 
teraction between the hidden and visible sectors is extremely 
weak, pi and p2 are identical.) Then, for any operator O in 
our visible sector, trpO = (trpiO + trp20)/2. This means 
that, as long as the variables measurable for us are concerned, 
the state can be understood as a statistical ensemble composed 
of pi and p2- Therefore what we actually observe is the ex- 
pectation value for either pi or p2- Therefore, irrespective of 
whether the isocurvature perturbation is in the visible sector 
or in the hidden sector, we can insert a projection operator V 
to take into account the influence of decoherence. In the suc- 
ceeding section, we discuss the regularity of the "observed" 
n-point function {V 0{tf, Xi)0{tf, X2) ■ ■ ■ 0{tf,Xn)). 



III. PROOF OF IR REGULARITY 

In this section, we prove the IR regularity of the n-point 
function Eq. ( |2.42| i. In this paper, we discuss the evolution 
of perturbation during a finite period of inflation. First we 
describe the way of the quantization in i jlll Al Before starting 
the detailed discussion, in ijlllBI we briefly explain the basic 
idea of the proof of IR regularity. In this subsection, we clarify 
the difference between the regularization in multi-field models 
and that in single field models. After that, in mil CI we adapt 
the basis transformation. In the new basis, it becomes easier to 
understand the regularization in the multi-field models. Based 
on these preparations, in i jlllDI we show that IR suppression 
due to the projection operator regularizes the IR divergence 
when the initial conditions are set at a finite past. 

In this section, we discuss the IR regularity after we remove 
the influence of the unobservable quantities. For the technical 
reason, it is better to avoid treating the divergent quantities di- 
rectly. Therefore, first, we assume that the total volume of the 
universe Vc = is finite. Then, the normal modes take dis- 
crete spectrum, and as a result IR divergence is concentrated 
on the spatially homogeneous mode with p = 0, as long as 
Lc is kept finite. Even with a finite volume, the quantum fluc- 
tuation of the homogeneous mode with p = in adiabatic 
vacuum is still divergent in contrast to the other IR modes. In 
Sec. nil Al we introduce a parameter Sa that measures the de- 
viation from the adiabatic one for the p = mode. At the end 
of calculations, we take the limit Vc ^ 00 and Sa 0. 

A. Quantization 

In the previous section, we described how we can expand 
the Heisenberg field 0^ in terms of (^{(a;). The interac- 
tion picture field ^{{x) ^ satisfies the equation of motion 
^iJ'Pi (x) — 0. Using a set of mode function {(l>^a,p{^) = 
p(t)e*P '^} which satisfies 

= e-p-/:Va,p 

(3.1) 

we expand <fi{{x) as 

^ii-) = ^EE |%^^"-^«-P+h-c.| , (3.2) 

where the index a — 1 • • • £> is the label of the orthonormal 
basis. Making use of the Gram-Schmidt orthogonalization, 
the mode functions (/)^ p = e'^'^u^ p/MpiVc^^ are orthonor- 
malized such that 

{4'a,p,4>P,p') = VcSapSpp' , (3.3) 



^ The leading terni of the Heisenberg picture field (fiKx) agrees with the 
interaction picture field ipf^_^f.(x). Thus, we denote the interaction picture 
field as ip^{x). 
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is satisfied, where the Klein-Gordon inner product is defined 
by 

(,/., V) = -m' / {(l^'dari - {da<l}') md^" ■ (3.4) 

In Eq. i3.3i the factor Vc is necessary in order that the same 
functional form of the mode functions satisfies the natural or- 
thonormal conditions in the continuum limit, Vc ^ oo. (See 
Appendix |A1 ) Using the normalization conditions ( 13.31 ), we 
find that the creation and annihilation operators aj^ p and a^ p 
satisfy the following commutation relations, 

[tta.p, CLp p,] = SapSpp' . (3.5) 

The initial vacuum state |0)a is annihilated by the operation 
of any annihilation operator: 

O'a,p\0)a = 0, for and . 

The mode function ^{t) is normalized by 

ui,p{t)u*ip,p{t) - uLp{t)u*ip,p{t) = 5a,l3 . (3.6) 

In the long wavelength limit, we obtain two real independent 
growing and decaying solutions as 

gi,p{t)=gi{t)[l + 0{{p/aHf)] , 

dL.pit) = diit) [1 + 0{{p/aHr)] , (3.7) 

and g'^ (t) and d^^ (t) satisfy the normalization condition 

gi{t)dif3{t) ~ gi{t)di0{t) = a-^it) S^^p . (3.8) 

{gen dp) are the time dependent solutions. In massless de 
Sitter approximation, they are given by = (5^ and 

d^^{t) ~ l/{&Ha?)5^^, respectively. Combining these two 
solutions, we can construct a mode function as 

uLk(i) = ^5 a,feW + *<(fc)d^„,fe W, (3.9) 

with an arbitrary parameter Ca{k). The squared amplitude of 
fc(^) gives the amplitude of the primordial perturbations. 
It is common to set the initial state to the adiabatic vacuum 
which is a natural state in the inflationary universe. At the 
horizon crossing, where k « aH, the growing and decaying 
solutions should contribute to the positive frequency function 
""q fc(0 '^o the same order unless the initial quantum state is 
very different from the adiabatic vacuum one. Assuming that 
the time variations of and a^Hd^^ f. are not very fast after 
the horizon crossing time, i.e. \g^^ k\/\d^a k\ ~ Ha^'^^" with 
(5q ^ 1, this requirement determines the order of magnitude 
of Ca{k) as 

Thanks to the local gauge conditions, as in the single field 
case discussed in our previous paper Ir42ii1 . we can prove the 



regularity of the IR fluctuations initially in the adiabatic direc- 
tion, i.e. the tangential direction to the background trajectory. 
Looking at Eq. ( 13.1b . we find that (p^ / p = d(j)^ /dp satisfies the 
mode equation for the homogeneous mode u^^ q. We choose 
one of the bases m( j, w ^./ci{k) so as to approach d(f)^ /dp 
in the homogeneous limit fc — > 0. Then, as we will show in 
i jlllBI the modes with a = 1 no longer cause IR divergences. 
We give the other modes ^ (a — 2, ■ ■ ■ D) so as to be or- 
thogonal to each other 

As we are considering the universe in a finite box, wave 
numbers k are discrete. Hence, unless we take the infinite 
volume limit, Vc oo, the divergence is concentrated on the 
spatially homogeneous mode with fc = in the above expres- 
sion for the mode functions. To deal with this divergence in 
the fc = mode, we regularize Cq(0), introducing a small 
parameter Sq,, as 

C„(0) = Sc^/Vc^ . 

Then, we obtain 

<,oi^) = — 5a,o(i) d'a.oit) , (3.11) 

Sa Vc" 

After we define appropriate observables, we take the limit 
Sa ^ and Vc ^ oo. 

Giving the Wightman function G^j{x, x') in Eq. ( 12.25b as 
G^j{x,x') —a (0|(^/^i(a;)(^j_i(a;')|0)a, we can expand the 
retarded Green function in terms of mode functions u^^ as 

Gr^x, x') = -i6{t-t') e*(^-^')i?S,fc(i,^'), 

" k 

(3.12) 

where 

D 

R\kit,t') ^ T.^<,ki*)''Jo.At')-c.c. }. (3.13) 

Q = l 

Then, using the expressions in Eq. ( 13. 9b . we find that i? j ^, is 
regular even in the IR limit k 0. 

B. IR vanishing smooth function 

In this paper we do not consider the secular growth of the 
amplitude of perturbation due to the integration for a long pe- 
riod of time. Namely, we consider the case that ti is set at a 
finite past from tf. Deferring the detailed explanation to our 
succeeding paper, we give a brief comment on the regulariza- 
tion of the secular growth in multi-field model in SeclIVI In 
this paper we concentrate on the IR divergences originating 
from the momentum integration. 

The first part of our proof of IR regularity in multi-field 
model goes in parallel with the single field case f42l]. In the 
single field model, the proof of IR regularity was quite sim- 
ple if we do not care about long time integration. However, 
multi-field extension turns out to be non-trivial even for this 
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restricted case. To keep the simplicity of notation, we sup- 
press the field indices / and the labels of modes a for a mo- 
ment. As (f{x) is composed of (n = 1, 2, 3, • • • ), we make 
use of the mathematical induction to show the regularity of all 
(fn- 'f7i{x) is, by definition, n-th order in the interaction pic- 
ture field ifi. Formally, we define C[(^„](a;;pj, • • • ,p„) by 
expanding (fin{x) as 



C(°)[(^„](x;pi,- • • ,p„) 



_j=i Vc p^^a Pj 

j = l Vc'' p^^oPj 

xCW[(^„](x;Pi,-- - ,p„_i) + --- ,(3.14) 



3/2 



where we have suppressed the terms containing creation op- 
erators. C^^-' represents the coefficient of the term which con- 
tains j-th order product of 0-mode operators oq. We also sup- 
press this superscript (j), for simplicity. The above expression 
is the result that we obtain after conducting all the integrations 
over the intermediate vertexes. The momenta {Pj} in the ar- 
gument of C[i^„] are those associated with the right most ends 
of the corresponding tree-shaped graph. 

A key ingredient of the first part of our proof is to show that 
C[ipn] {x;Pi, - ■ ■ ,Pn) has the following properties, 

• It is a smooth function with respect to x for = 
\Pj\ < a{t)A, where A is an UV momentum cutoff 



scale. 

• It vanishes when the long wavelength limit pj 
taken for any momentum in its arguments. 



is 



If C[(p„] satisfies the properties mentioned above, (then 
we say C[i^n] is an IR vanishing smooth function 
(IRVSF)), one can easily show that 7i-point functions 
{(f{tf,Xi)--- if{tf,Xn)) are free from IR divergences. When 
we take the expectation value of the product of cpj {j < n) in 
the form of Eq. (|3.14| i, we consider all the possible ways of 
pairing a„. with aj, . Then, each pair of a„. and a], is re- 
placed with 6p _p . One of the momentum integrations over 
Pj and Pj is performed to obtain an expression in the form 



-^C[^„J(xi; • • • ,Pj, • • •)C[(^„J(x2; • • • ,Pj,---). 



in the continuum limit. Here we note that Vc ' J2j=i '^Pj 
should be replaced with {2tt)^'^ J (Ppj Up in the continuous 
limit. (See Appendix lAl) The resulting momentum integra- 
tion does not have IR divergences owing to the second prop- 
erty of C[(fn], i.e. limp^o C[(^„] (x; • • • , p, • • • ) = 0. 

Before we start the mathematical induction, let us note the 
following properties of IRVSFs: 

Lemma If Ci{x;{pj}) and C2{x\{qj}) are IRVSFs and 
there is no overlap between the list of momenta 
{Pj} and {q^}, then V aCi{x; {pj}), xCi{x;{pj}), 

C,{x;{p,}), A-iCi(x;{p^.}), WtCi{x;{p^}), 



J dtCi{x;{pj}), and Ci{x;{pj}) x C2(x;{qj}) are 
all IRVSFs. 

Now, let us prove that C[(p„] is IRVSF by induction if 
C[(fi] is so. The n-th order perturbation is obtained by 



= Wt dt' d^x'a^t')GR{x,x')Wt'{x')rn{x') . 



(3.15) 

Wt' (a;')r„(.T') is constructed from lower order perturbations 
SNj, Xj, fj, B±j and (^j with j < n using the operations 
listed in the above Lemma. Furthermore, from Eqs. ( 12.281 1, 
(|230] | and ( |236] |, we find that SNj, Xj and fj are all 
constructed from ifi with I < j hy the operations listed 
there, too. Hence, C[M^t'r„], the expansion coefficient of 
Wt'{x')rn{x') analogous to C[(p„] in Eq. ( 13.141 ). is also an 
IRVSF. Since the Fourier mode of the retarded Green func- 
tion described by Eq. ( 13.131 1 is regular in the IR limit, its 
Fourier transform Gii{x,x') should be regular, too. (Reg- 
ularity in UV is assumed to be guaranteed by an appro- 
priate UV renormalization.) Since the integration volume 
of x' is finite, the integral of a product of regular func- 
tions Jdh'a^{t')GR{x,x')Wt'ix')Tnix') should be finite, 

and hence it is IRVSF. Since the operation Wt preserves the 
properties of IRVSF, — TytGfl(Wt'r„) is also foundtobe 
IRVSF 

Now our concern is whether the first step of induction is 
true or not. Namely, we examine if C[(^i] is IRVSF or not. 
Utilizing the residual gauge degrees of freedom, we fix the lo- 
cal average of the adiabatic mode A = ejCp^ . Then, IR modes 
in this direction are controlled to be free from divergences, 
but the modes in the other directions are not. The interaction 
picture field appears only in the projected form Wjip'l, which 
can be expanded by the mode function u^^ j, as 



w'jf>iix) 



where 



Vc 



Wt-r 



^ipx 



-{h.c.}, 



Wt 



p I 

e ej 



pi 



d^x e'P '' Wt{x) 



(3.16) 



(3.17) 



and we note that Wt^o = J d\ Wt [x) — L^.To make it easy 
to take the limit Vc oo, we define the Fourier mode of the 
window function in a different manner from those of fluctua- 
tions. (See Appendix IaI) Hence, we have the coefficient for 



C'^^\w'jCpi]{x,p)^ 



ip-x cl 
3 ^ 



Wt. 



Wt 



-p I 
— e ej 



it) 



Mpi 
(3.18) 

We have chosen the adiabatic mode (a = 1) so as to be 
tangential to the background trajectory in p ^ limit, i.e., 
g\{t) (X e^. Then, multiplying p^u-^j^ p(t) = by the 

factor [e'P-^ - Wt,-p/Wo], c["\w^j^(]{x,p) vanishes in 
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this limit. Therefore Ca\Wjip{]{x,p) vanishes in the limit 

p 0. Thus, we find that da'^[W^j(p(]{x,p) is an IRVSF. 
However, the factor [e^P '^S^j — Wt^-p/WtfiC^ej] does not 
suppress isocurvature fluctuations , which is pointing or- 
thogonal direction to the background trajectory. When one of 
the basis with a ^ 1 has non-negative value of 6a, the IR con- 
tribution of such a mode diverges and Ca'' [Wj(f(]{x,p) is 
not IRVSF. In this case n-point functions of 0{x) — {A, S^} 
actually diverge. 

The case with p = goes in a similar manner. The coeffi- 
cient for Ua.o of (^{ is given by 



gjoit) 

SaMpi 



(3.19) 



where we take the limit Vc ^ oo using Eq. ( 13. lit . This ex- 
pression also vanishes for the adiabatic perturbation, but it 
does not for the isocurvature perturbation. Then, the contri- 
bution from the isocurvature perturbation diverges when the 
limit Sq, —> is taken. 

The above divergences arise only in the multi-field model. 
In the rest of this section, we discuss the regularization of this 
divergence. Even when 6a is negative, the following discus- 
sion is still relevant. For (5g < there is no IR divergence, but 
IR contribution can be large if |(5a | <C 1. If one can show that 
n-point functions for 0{x) are free from the IR divergence for 
6a > 0, it also implies the absence of enhanced IR contribu- 
tion for 6a < 0- 



C. Squeezed wave packet 

For the later use, we transform the mode functions 
{u^a '^o another ones suitable for discussing the effect 

of inserting projection operator Transformation proceeds in 
two steps. Deferring the detailed explanations to AppendixiBl 
here we just give a brief sketch of the transformations. At the 
first step the new basis mode functions {v^a p} for P 7^ are 
arranged so that the leading term in the long wavelength limit 
is cancelled. While, the IR divergent contributions are local- 
ized to a single mode with p = 0. At the second step of trans- 
formation, without changing the mode function for p 7^ 0, we 
introduce another mode function for p ~ mode which natu- 
rally defines wave packets with a finite width even in the limit 
Sa ^ and Vc ^ 00. In this limit, ip{ can be expanded as 



where 



a-) = E ■ 



a,0 



aa,0 + 



+ (h.c.) , (3.20) 



t I (fp Wp I 



Hf Jp^o (2ny Wo 



v' (x) ^ u' me^P---^^^^ 



(3.21) 
(3.22) 



Now, we introduce the coefficients C^a} ["^ril analogous to 
C|^j ["^rJ foi' '^he creation and annihilation operators aj^ ^ and 
cia.p- Then, in the continuum limit, if^ is expanded as 



k=0 



.1=1 

^(fc) r 



n 



(fpi 



'■aj.pj 



(3.23) 



Since the induction with respect to n proceeds as before, we 



can say that cj^^j [^fj 



is IRVSF if the coefficients of the first 



order variables Wjipi{x) are IRVSFs. 

In the same way as Eq. ( 13.181 ), the coefficient for aa.p is 
obtained as 



Wt.-p I 
— — e/ 



M, 



pi 



(3.24) 



From Eq. ( 13.22b . p^^^v'^ ^ vanishes in the limit p ^ 0. The 
adiabatic component also vanishes owing to the projection 

W^j. Hence these coefficients are IRVSFs. While, the co- 
efficient for o is given by 

C^^\w'j(pi]{x) = [6'j 



it) 



M, 



pi 



(3.25) 



which is finite from Eq. (13.21b in contrast to the previous case 
in Eq. ( 13. 19b . Therefore it is also IRVSF. As we find that all 

the coefficients of W^jCp'l are IRVSFs or just a regular func- 
tion independent of p, all higher order coefficients cj^j [1^^] 
are proven to be IRVSFs by induction. 

We should emphasize that even if the coefficients c'"^^-^ [ipl^ 
are IRVSFs, this does not imply the regularity of n-point func- 
tions of O for the initial adiabatic vacuum state |0)a. In the 
following discussion we use an expression for the adiabatic 
vacuum state \Q)a expanded in terms of the coherent states 
1/3) a associated with ciafl- The coherent state satisfies 



aipl/ 



)a=(i&\ 

>a = 0. 



As shown in AppendixiBl the original vacuum state | )q can 
be expressed as 



|0)„ = 



D 

n 



dPaEaifia) 



where the coefficient Ea{Pa) approaches to 



EaiPa) 



(3.26) 



(3.27) 



in the limit of 0. The non-vanishing support for Ea 

extends to infinitely large \(3a\ in this limit. This is a con- 
sequence of the fact that the wave function in the adiabatic 
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vacuum is highly squeezed in the direction corresponding to 

Using Eq. ( 13.26b . we can expand the n-point function of the 
"observables" 0{tf, x) as 

{rditf,xi)d{tf,x2)---d{tf.x„) ) 

D r 



n 

a=2 



x^{(3\Vditf,x,)---ditf,x^)\j)l° 



(3.28) 



where, to clarify that we should sum up only the connected 
graphs, we have added the suffix "conn". 

One remark is in order in computing the expectation value 
of the product of {oq, p} and {aj^, p}. Basically pairs between 

aa,p and a^, p, are replaced with S{p — p')5a,a' except for 
the case with p = 0. After the replacement, only the op- 
erators {aqo} and {ajj, q} are left on the right hand side of 
Eq. (13.281 1. The expectation value of the product of these op- 
erators becomes the summation of the commutators and their 
normal ordered products. Since the operators a^.n and a)^ q 
are not annihilated by the coherent state, the expectation value 
of the normal ordered products composed of these zero-mode 
operators are non-vanishing. The annihilation operators o 
in the normal ordered parts acting on the coherent state |7)a 
produces the factor 7^, while the creation operator q act- 
ing on a{P\ produces In other words, a^ o and a\ q are 
replaced either with a commutator by making a pair or with 
7q. and respectively. These two exclusive possibilities can 
be concisely expressed by the replacement 



la 



where o^'q and a|,'|/ satisfy the same commutation relation as 
(lafl and q, and they annihilate the coherent state |7)a and 
a{/3\, respectively. 

Now it will be obvious that n-point function evaluated for 
the coherent states, 

^{f3\d{tf,x,)---d{tf,x„) |7)r"- 

is finite. To show its regularity, the insertion of the projection 
operator V is unnecessary. We were focusing on C'j^^'j [1^^], 
which is the coefficients of the annihilation operators, but all 
the other coefficients of the mixture of {da.p} and {a^, p} are 
shown to be IRVSFs in the same manner. The effect of coher- 
ent state is taken care by the replacements (|3.29t . Therefore, 
the expectation value for fixed values of and 7a is regular. 



D. Role of projection operator 

Without the projection operator V, the n-point functions of 
0{x) for the adiabatic vacuum diverge, although the expecta- 
tion values for the coherent states were proven to be finite in 
the preceding subsection. The divergences for the adiabatic 



vacuum appear in the integration over f3a and 7q,. The limit- 
ing behavior of Ea{/3a) for Sa ^ given in Eq. ( |3.27t tells 
that this factor does not restrict the effective range of these in- 
tegrations in the limit 0. Therefore, integrating over /3a 
and 7q, without the insertion of V, the n-point function for the 
adiabatic vacuum diverges. This result is as expected since 
the basis transformation does not change the final result for 
the n-point function. 

In order to remedy these divergences, we need the insertion 
of the projection operator. The projection operator V takes 
care of the effect of quantum decoherence, removing the con- 
tamination from the other parallel worlds. We will see that the 
insertion of V makes the effective range of integration for /3a 
and ja finite. 

In the same way as 0{tf,x), we expand V, which is the 
functional of {S"{tf)}, in terms of {da.p} and {aj, p}. Ex- 
panding S°'{tf) as 



s''{tf)^sntf) + lsntf) 



(3.30) 



we focus on the coefficient of da,o + dl^Qin the leading term 
5f (t/) = eJ{Wtfif>{ + We decompose the terms 

which contain g and at g into two pieces proportional to 
{dafi + d\ q) and {da.o — at Then, we can see that the 
former does not contain contribution from S^j^ as follows. 
S^j^ is determined by solving Eq. ( 12.361 ), which is sourced by 
ejWt^pi on the right hand side. The coefficient of dafi + at q 
in Wt(p{ is given by 



(3.29) + vt,) = Hfgi, + f 



d^p 
(c.c.) ) .(3.31) 



WpWt,-p J 
WoWt.o 



Using the relation Wt,-p — W^p, which is derived from the 
reality condition of Wt{x), the second term on the right hand 
side vanishes. Contraction with e/ erases the first term, too. 
Since the source term of the equation for B^^^ does not contain 
{dafi + at q), B^j^^ does not, either. 

Using the replacement ( 13.291 ), we divide (ag o + at p) into 
(Pa + la) and (flg + a^.o^)- Separating the coefficient of 
(/3a + 7a), 5f (t/) is expressed as 



5f(t;) = Al(/3^>7'') + 55f(t/), 



where 



(3.32) 



(3.33) 



and summation over repeated index (3 is understood. 

Inserting the above expression for S°'{tf) into V given in 
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( |2.39l l, the observable n-point function ( 13.281) is recast into 

{V0{tf,xi)0{tf,X2)---0{tf 



D r 



n 



exp 



oo J — oo 



I exp 



(3.34) 



where 



Owing to the first exponential factor, the contribution from 
the integration region with | A1(7'' + ^ cr is exponen- 
tially suppressed. Since the inner product between the coher- 
ent states gives 



a(/3|7)a = flexp 



(see Eq. ( |B13l l), the contribution from the region with — 
7" I 3> lis also exponentially suppressed. The directions of 
these two suppression are orthogonal, and hence the effective 
integration area is restricted to a finite region 



< 



(3.35) 



where A and a are the typical amplitudes of the eigenvalues 
of A^ and the square of the eigenvalues of C^, respectively. 
Our discussion up to here ensures the finiteness of the effec- 
tive range of the Gaussian integrations over /3" and 7". This 
proves the IR regularity of the n-point function of the local 
fluctuation {0{x)} with the projection at each order of loop 
expansion even if the initial state is set to the adiabatic vacuum 
state. 

Now the remaining task is to examine if the perturbative 
expansion is still reliable after all the changes that we made. 
For a sufficiently wide projection operator V , the expected 
amplitude of \(p^\ is dominated by the contribution from /3" 
and 7", which is 0{(t). Hence, the validity of the perturbative 
expansion requires 



tJ < 1. 



(3.36) 



The next question is whether one can safely expand the second 
exponential factor in (|3.34t . 



coming from V. Using the conditions ( 13.351 ). we find that this 
expansion converges if 

is satisfied. On the other hand, the amplitude of is esti- 
mated to be given by the linear order contribution A. There- 
fore the necessary condition is 



A < a. 



(3.37) 



Since A ^ 1, we can choose a such that satisfies (13.36b and 
(13.37b simultaneously. 



IV. CONCLUSION 

In the present paper, we have proposed one solution to the 
IR divergence problem in multi-field models. We discuss n- 
point functions for the local perturbations variables, 0{x), de- 
fined by the deviations from the local average values with an 
additional gauge condition that fixes one of the residual gauge 
degrees of freedom remaining in the usual flat gauge. Even 
if we consider these local perturbative variables, when plural 
fields have a scale invariant or red spectrum, we encounter IR 
divergences. This is because the effects of quantum decoher- 
ence are not taken into account yet; the n-point function for 
0{x) is affected by the contaminations from other uncorre- 
cted worlds. To remove the contaminations, we have inserted 
the projection operator V which projects the final quantum 
state to a wave packet with a sharp peak in the space of the 
local average values of the fields, {WtfSj}. Here, we give 
an intuitive way to understand how the insertion of the pro- 
jection operator regularizes the IR corrections. When plural 
fields contribute to IR divergences, the wave function corre- 
sponding to the initial adiabatic vacuum is highly squeezed in 
the corresponding directions in the space of {Wt^Si}. How- 
ever, only a part of the squeezed wave function does contribute 
to a decohered wave packet, which represents our world. In- 
troducing the projection operator, we have taken into account 
the restriction of the wave function to the non-vanishing sup- 
port of the decohered wave packet. This restriction is recast 
into the exponential factor in Eq. ( 13. 34b , by which the non- 
vanishing support of 7) becomes a finite region. This as- 
sures the regularity of the observable n-point functions. 

The "observable" n-point function (12.42b depends on the 
parameter a that we introduced to incorporate the decoher- 
ence effect without discussing the detailed process. It also 
depends on the size of the observable region. These depen- 
dencies may disappear when we compute the actual observ- 
ables, taking into account the decoherence process appropri- 
ately. Definitely, to predict the accurate value of observables, 
further study is necessary. We leave this issue for our future 
work. 

Here, fixing the temporal coordinate on each vertex, we 
showed the regularity of the integration over the spatial co- 
ordinates. Hence, we cannot deny the possibility that the tem- 
poral integral makes the n-point function diverge when we 
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send the initial time ti to a distant past. In single-field mod- 
els, we showed, by considering the local fluctuations, that 
we can suppress the secular growth which appears from the 
temporal integration, unless very higher order perturbation is 
considered. We expect that in the multi-field model, a simi- 
lar suppression appears from the inserted projection operator. 
The reason is as follows. If the non-vanishing contributions 
from the distant past makes the amplitude of 0{tf, x) sig- 
nificantly large, they must increase also the amplitude of the 
local average {WtfSj} to a large extent. However, the deco- 
hered wave packet we pick up has the bounded average val- 
ues of {WtfSj}. Therefore, we expect that as long as we 
compute the observable n-point function evaluated for a de- 
cohered wave packet, the contributions from the distant past 
are suppressed and the temporal integration converges. We 
also defer the examination of our optimistic expectation to a 
future work. 

In general, the IR divergences originate from massless 
(or quasi-massless) fields with non-linear interactions in the 
(quasi-) de Sitter universe ||llilll59, 60, 611111. In flSlTt], 
using the prescription with an IR cutoff, the logarithmic am- 
plification is discussed for a massless test scalar field with a 
quadratic interaction term. We can discuss the regularization 
of the IR corrections from a test field in a similar manner to the 
entropy perturbations discussed in this paper. In [59, 60], the 
effects of IR gravitons which grow logarithmically are argued 
to screen the cosmological constant. However, the Hubble pa- 
rameter defined in |59, 60] is gauge dependent and in other 
definition we do not encounter the screening 16 1|] . Thus, this 
problem is still controversial ll6ll l62ll . We would like to apply 
our prescription to these issues in our future work, too. 
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APPENDIX A: CORRESPONDENCE IN THE CONTINUUM 
LIMIT 



We use two different notations for the Fourier components be- 
tween the perturbation variable like f''{x) and the window 
function Wt^X^)- Because of that, these two classes of quan- 
tities should be treated differently when we take the contin- 
uum limit. Hereafter, we use variables Q and W to represent 
the variables of the first and the second classes, respectively. 
We add the index (d) on the discrete Fourier components in 
the finite volume, while the index (c) on the continuous ones 
in the infinite volume. If not necessary, we abbreviate the suf- 
fixes, / and a. 



1. Quantized variables 

When we consider the first class of variables Q{x), which 
is to be quantized like (p{x), the corresponding mode func- 
tions {qk} play the more important role than the Fourier com- 
ponents of Q{x) themselves. Therefore we adopt a conven- 
tion so that Qk remains unchanged in the continuum limit 
Vc oo. 

When we quantize Q{x), we expand the Fourier mode Qk 
in terms of the creation and annihilation operator hke 



Qk = (Stt) 2 {akQk + alql} 



(Al) 



We require that the mode functions qk remain unchanged in 
the continuum limit. On the other hand, the commutation re- 
lation for the creation and annihilation operator 

■ (d) (d)tl _ c 

changes in the continuum limit to 

6^{k-k') . 

From the above commutation relation, a trivial relations 



(c) (c)t 



1 = 



El' 



id) „(d)t 



k ' "'k' 



1 = d'k' 



follow. Since the wave number k is discrete like ki = Akji = 
2TTji/Lc with ji being integer, J^k' '■° replaced with 
{27r)~^Vc J d^k in the continuum limit. This requires the cor- 
respondence like 



Ad) 



V~ 



^k ' 



where Vc = Vc/(27r)'^ and hence we find 



Vc'"^Q^k^- 



(A2) 



In this paper, to tame the divergence in the IR limit, we 
begin with the model with a finite volume of the universe 
Vc = L^c- After we define appropriate observables, we take the 
infinite volume limit ^ oo. At this step the discrete label 
of the co-moving wave number changes to the continuum one. 



To realize the above correspondence, we define the Fourier 



components Q^^^ as 

^k — - 1 



d^xQ{x)e 



'ik-x 



(A3) 
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Then the definition of the Fourier components in the contin- 
uum limit becomes a normal one: 



cPx Q{x)e 



The inverse transform is given by 



Q{x) 



1 



?(<!) ik-x 
k ' 



and its continuum limit consistently recovers 
1 



Qix) 



1 



(d) ik-x 
k ^ 



(FkQ 



(c) ik-x 



(A4) 



(A5) 



(A6) 



(A7) 



(27r)3 ' - --"k 

In ^BT] we treat the fc = mode separately. Even in the 
limit Vc — » oo, this mode is treated as a discrete spectrum. 
Hence, the corresponding commutators keep the normaliza- 
tion in terms of Kronecker 5, and we have 



" Q.O' " Q.O 



1. 



This means that the correspondence of the creation and anni- 
hilation operators should be like a'^^'^Q ^ ^''^V Taking into 
account the exception for the fc = mode, Eq. ( |A7t is to be 
modified to 



Qa{x) 



a,0 



1 



(A8) 



(27r)l (27r)3 

where we have denoted the contribution from the basis a. 
Comparing this with the discrete case ( IA6b . we find that we 
should have 

Correspondingly, we have 



9a0 



I Ad) 

1 a.O- 



Finally, we mention the normalization conditions of mode 
functions. The normalization conditions for q^^^ given in 
Eq. (|3.3l l are understood as 

M) Jd)\ = V,Skk'= [ d^xe'^"-"")-^. (A9) 



Hence, in the limit Vc — > oo the right hand side is to be un- 
derstood as (27r)3(53(fe — k'). Therefore the normalization 
conditions of mode functions in the continuum limit should 
be 



(AlO) 



To summarize, the correspondence between the discrete 
Fourier components and the Fourier components in the con- 
tinuum limit is given by 



Vc^Q^k\ 



(c) 

Ik 



(d) 

<l k ' 



Vc^a'j^\ (All) 



for all modes k with a ^ a and k Q with a ~ a and 



Q 



Vc 



Q,0' y a.O '^c 1 a.O' " Q.O 



Ad) 

^ Q,0- 

(A12) 



2. Un-quantized variables 

The quantities of the second class, W{x), are not supposed 
to be quantized. In this case it is more convenient to con- 
sider the Fourier components that remain unchanged in 
the continuum limit: 



id) 



For this purpose, we simply define the Fourier components 
in a usual manner by 



"^k^ = / d^xWix)e 



— ik-: 



Then the inverse transform is given by 



k 



id)jk-. 
k ^ 



(A13) 



(A14) 



and its continuum limit Vc oo becomes 

k 

1, ^ 

k 



1 



(2^) 



d^kwi'^^e'''' 



(A15) 



as is expected. 



APPENDIX B: BOGOLIUBOV TRANSFORMATION 



1. Another set of mode functions 



In Sec. IIII CI we transformed the mode functions to more 
suitable ones to understand the role of the projection operator 
In this subsection, we give the detailed explanations for the 
two transformations. On the first transformation, the new ba- 
sis mode functions {tj^^ p} are constructed so that the leading 
amplitude in the long wave-length limit is canceled forp ^ 0. 
While, the IR divergent contributions are localized to a single 
mode with p — 0. Hence, when we expand (^j^ in terms of the 
creation and annihilation operators associated with the new 
set of mode functions, the expansion coefficients C[(^^J are 
IRVSFs except for the case with pj = 0. Such transformation 
can be achieved by taking 



4(0) 



p#0 



clip) 



(Bl) 



where we denote Wtf,p by Wp. The normalization constant 
Ma is chosen as 



W.p Ca{Q) 
Wo Ca{p) 



= l + 0{Sa) 



(B3) 
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where the second equality immediately follows from the ob- 
servation that only the term with p = remains in the limit 

Sa -> 0. 

The orthonormal relation for {v^^ ^} is given by 



{v Q,o,f/3,o) 



VcSafS, 

0, 



|M^oPc„(p)c*(p') j 
K'5a/3{<5p,p' + 0(4/V;)}, 

{V a,0,V*p .p) = {v a,p,vl ,p,) . (B4) 



The mode functions v^^ are not mutually orthogonal be- 
fore we take the limit ^ in Eq. ( IB4b . After taking the 
limit Sa — + 0, however, {v^^ p) becomes a set of orthonor- 
mal bases. We denote the creation and annihilation operators 
associated with {v^^ pi p ^^'^ ba.p, respectively. 

After this transformation, the IR divergent contribution is 
confined to v'^ o(^)- Indeed, the problematic term in u^^ p 
which scales as p-(-^+^'>)/^ in the IR limit is cancelled by 
the second term in Eq. (IB2b . Then, even when Sa > 0, as 
long as Sa < 1, we find that pi 



a.pT^O 



vanishes in the limit 



p 0. This is the necessary and sufficient condition for the 
coefficients of the creation and annihilation operators to be 
IRVSFs. In contrast, v^^ q has a diverging amplitude in the 
limit Sa and Vc ^ oo. Notice that this transformation 
does not mix the positive frequency modes with the negative 
frequency ones. Therefore in the limit Sa ^ the vacuum 
annihilated by ba.p is identical to the vacuum annihilated by 

Now we move on to the second step of transformation. Here 
we introduce another mode function for p = mode which 
naturally defines wave packets with a finite width even in the 
limit — > and 14 ^ oo. In Fig.IT] we depicted the changes 
of the wave packets under these two Bogoliubov transforma- 
tions. For brevity, we have suppressed the indices / and a. 
We denote a new set of mode functions by {v^^ p}, which is 
defined by 

^^a,oi^) = coshro, w-^„_Q(a;) - sinhra 0(2;), (B5) 
v'aA^) = ^'a,ki^), (B6) 
with the squeeze parameter such that 

1 



(B7) 



where Aa is a real constant. The mode functions for p 7^ 
are unchanged: 



for p ^ 0. 



(B8) 



We denote the creation and annihilation operators associated 
with {v^a,p} by aj,p and dap- 
In the limit s„ 0, we obtain 



pk 













>■ 


< 


>■ 

, -3/2 



Uk 



1 









>- 






^ 


T2 


J 


>■ 

~ 1/ s 

i 



Vo 



^ Vo 



~A 

FIG. 1: Two Bogoliubov transformations 



p/O 



(B9) 



The first and the second terms in Eq. 



originate from the 



homogeneous mode. 



and its complex conjugate. The 



terms that diverge in the limit Sq, — > are common in u^^ q 
and u^^ g*, which cancel with each other in the definition of 
v^a ill 1^1- < IB5b . The third term in Eq. ( IB9I ) originates from 
the inhomogeneous modes, u^ap^o- The contribution from 
u^^ p^o i^ ^1^° large in g in the IR limit, but it is not in 
The leading terms of d^^ ^ (i) in p ^ limit behaves 



as 



Hence, the momentum integral J d^p Wpd^a p(^)6*''^' which 
arises in the continuum limit, remains finite for x < 00. 

Here, we choose the unspecified constant Aa such that the 
mode functions Va take the minimum amplitude on aver- 
age. To do so, we compare the amplitude of the first term 
with that of the last term in Eq. ( |B9b . Near the final time 
tf the last term is negligibly small, but it is larger in the 
past. Assuming and H are almost constant, the relative 
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amplitude between these two terms is roughly estimated as 
\^a^ lo" 'i^pd'a.pl ~ ^^Mq' where we used d'^p w 
l/{6Ha^)5{^. Hence, we find that it is appropriate to choose 
Aa = H to minimize the amplitude of the mode functions. 
(We choose Ai = 1/si = oo so that the mode function for 
the adiabatic mode is unchanged: v{ q = v{ q.) 

Now, sending Vc oo, we take the continuum limit. We 
have summarized the correspondence between the discrete de- 
scription and the continuous one in Appendix lAl The overall 

factor Vc^ in Eq. ( IB9b is cancelled by Vc ^ in Eq.( l3.2b . In this 
limit, replaced with {2tt)^'^Vc J cPk. Since the expres- 

sion iB% has a factor 1/ Vc in front of the summation 
find that v^^ q is free from the explicit divergence in the limit 
— > oo. In the continuum limit, can be expanded as 



a,0 



aa,0 + 



+ (h.c.) , (BIO) 



where 



Jp^o (^T^)' Wo 



Wo Caiv) 



(BID 
(B12) 



2. Expansion by coherent states 

In this subsection, we show that the adiabatic vacuum state 
0)a can be expanded in terms of the coherent states |/3)a 



associated with a^.o- The coherent state \P)a is defined by 



a=2 



a=2 



(B13) 



where | )a is the vacuum state annihilated by aa,k- This new 
vacuum state | )a is related to | )a by 



I )a = n (cOShra)-^/2gitanh., (4, J= | g ^ (3^4^ 

a=2 

It is easy to verify that the left hand side is annihilated by 

ba.a — (cosh rg) ag o ~ (sinhra) flg Q. The coherent state 
satisfies 



ai.p\l3)a = 



a IK/ a J 



The original vacuum state \0)a can be expressed as a super- 
position of the new coherent states like 



^ r /'OO 
|0)a = n / dP^E^ilSo,) 

where the coefficient Eg, {Pa ) is given by 

EaiPa) = (2^sinhra)-i/2g-/3l(e--~i)- 



(B15) 



(sa->0). (B16) 



The formula dBlSb can be verified simply by performing 
Gaussian integral about {/3a}. 
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